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Systems with long-range interactions display a short-time relaxation towards 
Quasi-Stationary States (QSSs), whose lifetime increases with system size. 
With reference to the Hamiltonian Mean Field (HMF) model, we here review 
Lynden-Bell's theory of "violent relaxation" . The latter results in a maximum 
entropy scheme for a water-bag initial profile which predicts the presence of 
out-of-equilibrium phase transitions separating homogeneous (zero magneti- 
zation) from inhomogeneous (non-zero magnetization) QSSs. Two different 
parametric representations of the initial condition are analyzed and the fea- 
tures of the phase diagram are discussed. In both representations we find a 
second order and a first order line of phase transitions that merge at a tri- 
critical point. Particular attention is payed to the condition of existence and 
stability of the homogenous phase. 
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1. Introduction 

Hamiltonian systems arise in many branches of applied and fundamental 
physics and, in this respect, constitute a universal framework of extraordi- 
nary conceptual importance. Spectacular examples are undoubtely found 
in the astrophysical context. The process of hierarchical clustering via grav- 
itational instability, which gives birth to the galaxies,^ can in fact be cast 
in a Hamiltonian setting. Surprisingly enough, the galaxies that we observe 
have not yet relaxed to thermodynamic equilibrium and possibly corre- 
spond to intermediate Quasi-Stationary States (QSSs). The latter are in a 
long-lasting dynamical regime, whose lifetime diverges with the size of the 
system. The emergence of such states has been reported in several differ- 
ent domains, ranging from charged cold plasmas^ to Free Electron Lasers 
(FELs),'^ and long-range forces have been hypothesized to be intimately 
connected to those peculiar phenomena. 

Long-range interactions are such that the two-body interaction poten- 
tial decays at large distances with a power-law exponent which is smaller 
than the space dimension. The dynamical and thermodynamical proper- 
ties of physical systems subject to long-range couplings were poorly under- 
stood until a few years ago, and their study was essentially restricted to 
astrophysics (e.g., self-gravitating systems). Later, it was recognized that 
long-range systems exhibit universal, albeit unconventional, equilibrium 
and out-of-equilibrium features.'' Besides slow relaxation to equilibrium, 
these include ensemble inequivalence (negative specific heat, temperature 
jumps), violations of ergodicity and disconnection of the energy surface, 
subtleties in the relation of the fluid (i.e. continuum) picture and the par- 
ticle (granular) picture, new macroscopic quantum effects, etc.. 

While progress has been made in understanding such phenomena, an 
overall interpretative framework is, however, still lacking. In particular, 
even though the ubiquity of QSSs has been accepted as an important gen- 
eral concept in non-equilibrium statistical mechanics, different, contrasting, 
attempts to explain their emergence have catalysed a vigorous discussion 
in the literature."' 

To shed light onto this fascinating field, one can resort to toy models 
which have the merit of capturing basic physical modalities, while allowing 
for a dramatic reduction in complexity. This is the case of the so-called 
Hamiltonian Mean Field (HMF) model which describes the evolution of N 
rotators, coupled through an equal strength, attractive or repulsive, cosine 
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interaction.^ The Hamiltonian, in the attractive case, reads 

JV JV 

^=2E^'' + ^E[l-^«<^.-^0], (1) 

where 9j represents the orientation of the j-th rotator and pj stands for 
the conjugated momentum. To monitor the evolution of the system, it is 
customary to introduce the magnetization, an order parameter defined as 
M = |M| = |^mi|/A^, where m.i = (cos 0,, sin 6*,) is the magnetization 
vector. The HMF model shares many similarities with gravitational and 
charged sheet models^'^ and has been extensively studied^ as a paradig- 
matic representative of the broad class of systems with long-range interac- 
tions. The equilibrium solution is straightforwardly worked out^ and reveals 
the existence of a second-order phase transition at the critical energy den- 
sity [/c = 3/4: below this threshold value the Boltzmann-Gibbs equilibrium 
state is inhomogeneous (magnetized). 

In the following, we shall discuss the appearence of QSSs in the HMF 
setting and review a maximum entropy principle aimed at explaining the 
behaviour of out-of-equilibrium macroscopic observables. The proposed ap- 
proach is founded on the observation that in the continuum limit (for an 
infinite number of particles) the discrete HMF equations converge towards 
the Vlasov equation, which governs the evolution of the single-particle dis- 
tribution function (DF). Within this scenario, the QSSs correspond to sta- 
tistical equilibria of the continuous Vlasov model. As we shall see, the theory 
allows us to accurately predict out-of-equilibrium phase transitions separat- 
ing the homogeneous (non-magnetized) and inhomogeneous (magnetized) 
phases. Special attention is here devoted to characterizing analytically 
the basin of existence of the homogeneous zone. Concerning the structure 
of the phase diagram, a bridge between the two possible formal settings, 
respectively^" and,^^ is here established. 

The paper is organized as follows. In Section 2 we present the continuous 
Vlasov picture and discuss the maximum entropy scheme. The properties of 
the homogeneous solution are highlighted in Section 3, where conditions of 
existence are also derived. Section 4 is devoted to analyze the stability of the 
homogeneous phase. A detailed account of the phase diagram is provided in 
Sections 5 and 6, where the case of a "rectangular" and generic water-bag 
initial distribution are respectively considered. Finally, in Section 7 we sum 
up and draw our conclusions. 
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2. On the emergence of quasi-stationary states: Predictions 
from the Lynden-Bell theory within the Vlasov picture 

As previously mentioned, long-range systems can be trapped in long-lasting 
Quasi-Stationary-States (QSSs),^^ before relaxing to Boltzmann thermal 
equilibrium. The existence of QSSs was firstly recognized with reference 
to galactic and cosmological applications (see'' and references therein) and 
then, more recently, re-discovered in other fields, e.g. two-dimensional tur- 
bulence^'* and plasma-wave interactions.^ Interestingly, when performing 
the infinite size limit N oo before the infinite time limit, t — > oo, the 
system remains indefinitely confined in the QSSs.-'^^ For this reason, QSSs 
are expected to play a relevant role in systems composed by a large num- 
ber of particles subject to long-range couplings, where they are likely to 
constitute the solely experimentally accessible dynamical regimes. ^''^ 

QSSs are also found in the HMF model, as clearly testified in Fig. 1. 
Here, the magnetization is monitored as a function of time, for two different 
values of N. The larger the system the longer the intermediate phase where 
it remains confined before reaching the final equilibrium. In a recent series 
of papers, '^'*''"*^'*^ an approximate analytical theory based on the Vlasov 
equation has been proposed which stems from the seminal work of Lyndcn- 
Bell.*^ This is a fully predictive approach, justified from first principles, 
which captures most of the peculiar traits of the HMF out-of-equilibrium 
dynamics. The philosophy of the proposed approach, as well as the main 
predictions derived within this framework, are reviewed in the following. 

In the limit of iV ^ oo, the HMF system can be formally replaced by 
the following Vlasov equation 

%+P%~ (M^m sin^ - My[f] COS0) |^ = 0, (2) 

where f{9,p,t) is the one-body microscopic distribution function normal- 
ized such that A4[f] = J fdOdp — 1, and the two components of the complex 
magnetization are respectively given by 

AUf] = J f cos edOdp, (3) 

My[f] = J f sin edOdp. 

The mean field energy can be expressed as 

+ Ml 1 



U 



\j fp^d9dp-'^^^^±^ + \. (4) 
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Fig. 1. Magnetization M{t) as function of time t. In both cases, an initial "violent" re- 
lajcation toward the QSS regime is displayed. The time series relative to N = 1000 (thick 
full line) converges more rapidly to the Boltzmann equilibrium solution (dashed hori- 
zontal line). When the number of simulated particles is increased, N = 10000 (thin full 
line), the relaxation to equilibrium gets slower (the convergence towards the Boltzmann 
plateau is outside the frame of the figure). Simulations arc carried on for a rectangular 
water— bag initial distribution, see Eq. (13). 

Working in this setting, it can be then hypothesized that QSSs corre- 
spond to stationary equihbria of the Vlasov equation and hence resort to the 
pioneering Lynden-Bell's violent relaxation theory-'^^ . The latter was ini- 
tially devised to investigate the process of galaxy formation via gravitational 
instability and later on applied to the two-dimensional Eulcr equation. 
The main idea goes as follows. The Vlasov dynamics induces a progres- 
sive filamentation of the initial single particle distribution profile, i.e. the 
continuous counterpart of the discrete N-body distribution, which proceeds 
at smaller and smaller scales without reaching an equilibrium. Conversely, 
at a coarse grained level the process comes to an end, and the distribu- 
tion function f(0,p,t), averaged over a finite grid, eventually converges to 
an asymptotic form. The time evolution of a rectangular water-bag initial 
distribution is shown in Fig. 2 

Following Lynden-Bell, one can associate a mixing entropy to this pro- 
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Fig. 2. The process of phase mixing is here illustrated, showing four snapshots of the 
time evolution of an initial rectangular water-bag distribution. The final state (right 
bottom) is a QSS. 

cess and calculate the statistical equilibrium by maximizing the entropy, 
while imposing the conservation of the Vlasov dynamical invariants. More 
specifically, the above procedure implicitly requires that the system mixes 
well, in which case, assuming ergodicity (efficient mixing), the QSS pre- 
dicted by Lynden-Bell, fQss{(),P,t): is obtained by maximizing the mixing 
entropy. As a side remark, it is also worth emphasising that the prediction 
of the QSS depends on the details of the initial condition,^^ not only on 
the value of the mass A4 and energy U as for the Boltzmann statistical 
equilibrium state. This is due to the fact that the Vlasov equation admits 
an infinite number of invariants, the Casimirs or. equivalently, the moments 
■Mn = J f^dddp of the fine-grained distribution function. In the following, 
we shall consider a very simple initial condition where the distribution func- 
tion takes only two values /o and 0. In that case, the invariants reduce to 
Ai and U since the moments Ain>i can all be expressed in terms of A4 
and fo as Mn = JJ^dedp = / f^""' x fdOdp = J f^'^JdOdp = f'o'~'M. 
For the specific case at hand, the Lynden-Bell entropy is then explicitly 
constructed from the coarse-grained DF / and reads 
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S[f] = -Jdpd0 



Jo fo \ fa J \ /o 
We thus have to solve the optimization problem 



(5) 



max{5[/] I U[f]^U,M[f]^l}. (6) 
/ 

The maximization problem (6) is also a condition of formal nonlinear dy- 
namical stability with respect to the Vlasov equation, according to the 
refined stability criterion of Ellis et al?^ (see also Chavanis^"'^). Therefore, 
the maximization of S at fixed U and M guarantees (i) that the statis- 
tical equilibrium macrostatc is stable with respect to the perturbation on 
the microscopic scale (Lynden-Bcll thermodynamical stability) and (ii) that 
the coarse-grained DF / is stable for the Vlasov equation with respect to 
macroscopic perturbations (refined formal nonlinear dynamical stability). 
We again emphasize that it is only when the initial DF takes two values 
/o and that the Lynden-Bell entropy can be expressed in terms of the 
coarse-grained DF /, as in Eq. (5). In general, the Lynden-Bell entropy is 
a functional of the probability distribution of phase levels. From Eq. (5), 
we write the first order variations as 

5S - (35U - aSM ^ 0, (7) 

where the inverse temperature (3 = l/T and a are Lagrange multipliers 
associated with the conservation of energy and mass. Requiring that this 
entropy is stationary, one obtains the following distribution^"'^^ 

/qSs(6',P) = ^ _^ e'3(pV2-A/x[/QSs]cose-M^[/Qsslsine)+Q ■ 

As a general remark, it should be emphasized that the above distribution 
differs from the Boltzmann-Gibbs one because of the "fermionic" denom- 
inator, which in turn arises because of the form of the entropy. Morpho- 
logically, this distribution function is similar to the Fermi-Dirac statistics 
so that several analogies with the quantum mechanics setting are to be 
expected. Notice also that the magnetization is related to the distribution 
function by Eq. (3) and the problem hence amounts to solving an integro- 
differential system. In doing so, we have also to make sure that the critical 



''The momentum P = f fpdOdp is also a conserved quantity but since we look for 
solutions where the total momentum is zero, the corresponding Lagrange multiplier 7 
vanishes trivially^^ so that, for convenience, we can ignore this constraint right from the 
begining. 
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point corresponds to an entropy maximum, not to a minimum or a saddle 
point. Let us now insert expression (8) into the energy and normalization 
constraints and use the definition of magnetization (3). Further, defining 
A = e" and m = (cos 0, sin 0) yields: 



/o J- /d0/_i/2(Ae-''^ ™)=l, (9) 



/o ^ (^^ ) / d0I,/, (Ae-^^-) = U 



fo 



J d0cos0/_i/2 (Ae-'3M m) ^ 
J d0sin0/_i/2 (Ae-''^ ™ 



fo ^^J d0sin0/_i/2 (Ae-'^^^^ ™) = M, 
where we have defined the Fermi integrals 



Ut) = I YZ^dx. (10) 







We have the asymptotic behaviours for t ^ 0: 



(-lnt)"+i 

Ut)-- f^, (n>-l), (11) 

n + 1 



and for t +00 



In{t)^ ^^''^^\ {n>~l). (12) 

The magnetization in the QSS, Mqss = Af[/Qss], and the values of 
the multipliers are hence obtained by numerically solving the above cou- 
pled implicit equations. It should be stressed that multiple local maxima of 
the entropy arc in principle present when solving the variational problem, 
thus resulting in a rich zoology of phase transitions. This issue has been 
addressed in^"'^^ and more recently in,^^ to which the following discussion 
refers to. 

It is important to note that, in the two-levels approximation, the 
Lynden-Bell equilibrium state depends only on two control parameters 
{U, fo) ^- This is vahd for any initial condition with f{9,p, t = 0) € {0, fo}- 
This general case will be studied in Section 6 where we describe the phase 



''These parameters are related to those introduced in'^'^ by (7 = e/4 + 1/2, /3 = 2r}, 
fo = no/N = Ai/(27r), k = 2tt/N, x = AO, y = (2/7r)Ap and the functions F in^^ are 
related to the Fermi integrals by Fk{l/y) = 2''°+-'-'/-^y7(fc_i)/2(y)- 
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diagram in the (/o, U) plane. Now, many numerical simulations of the N- 
body system or of the Vlasov equation have been performed starting from 
a family of rectangular water-bag distributions. The latter correspond to 
assuming a constant value /o inside the phase-space domain D: 

D = {{6,p) e [-7r,7r] x [-cx),c3o] | |6i| < M, \p\ < Ap}, (13) 

where < A6 < tt and Ap > 0. The normalization condition results in /g = 
l/(4A6'Ap). Notice that, for this specific choice, the initial magnetization 
Mq and the energy density U can be expressed as functions of A9 and Ap 
as 

M..=^, ,14, 

For the case under scrutiny, < Mq < 1 and U > Umin{Mo) = (l-M^)/2. 
The variables {Mq, U) arc therefore used to specify the initial configuration 
and hereafter assumed to define the relevant parameters space. This par- 
ticular but important case will be studied specifically in Section 5 where 
we illustrate the phase diagram in the {Mq,U) plane for the rectangular 
water-bag initial condition. Before that, we analytically study the stability 
of the Lynden-BcU homogeneous phase: two important limits, namely the 
non degenerate and the completely degenerate ones, are considered. We also 
discuss the condition for the existence of a homogeneous, non-equilibrium 
phase. 

3. Properties of the homogeneous Lynden-Bell distribution 

If we consider spatially homogeneous configurations {Mqss ~ 0), the 
Lynden-Bcll distribution becomes 

hssip) = Yri^im^- 

Using Eqs. (9), the relation between the inverse temperature (3 and the 
energy U is given in parametric form by 

''-«'V^/-v,(A)^ (16) 

This defines the series of equilibria T{U) for fixed /o parametrized by A 
(sec Fig. 3 in^"). The stable part of the series of equilibria is the caloric 
curve. Note that the temperature T is a Lagrange multiplier associated 
with the conservation of energy in the variational problem (7) . It also has 
the interpretation of a kinetic temperature in the Fermi-Dirac distribution 
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(15). If we start from a water-bag initial condition, recalling that fo 



l/(4A0Ap) and (Ap)^ = 6[C/ — (1 — Mq )/2], we can express /o as a function 



of Mo and U by 

f2 ^ I (lY) 

° 48[(2t/- l)(A6i)2 +sin2 A6']' ^ ' 

where A9 is related to Mq by Eq. (14). Inserting this expression in Eqs. 
(16), we obtain after some algebra the caloric curve T{U) for fixed Mq 
parametrized by A: 

^ 1 _ sin^A^ 



2 ^/-i/2(A)^ „2(Ag)2' 



6 I1/2W 

' ^^/-i/2(A)^-2(A^)^-^). (18) 



sin^ A6I V 6 /-i/2(A) 
Eqs. (16) can be rewritten 

(t^-J)8vrVo^ = G(A)^/l^, (19) 

where G(A) is a universal function monotonically increasing with A (see 
Fig. 2 of^"). A solution of the above equation certainly exists provided: 

(U - i)87r2/o' > G(0). (20) 

To compute G(0) we use the asymptotic expansions (11) and (12) of the 
Fermi integrals. This yields G(0) = 1/12. Therefore, the homogeneous 
Lynden-Bell distribution with fixed fo exists only for:^*^ 

1 1 
96^ 

For the rectangular water-bag initial condition, using Eqs. (17) and (21), 
we here find that the homogeneous Lyndcn-Bcll distribution with fixed Mq 
exists only for: 

1 / sin^Aei 



U > C/™„(/o) ^ T^^TTTI + o- (21) 



f/>t/™»(Mo).-^^,-^ + lj. (22) 

This result can also be obtained from Eq. (18) by taking the limit A ^ 0. 

Let us now describe more precisely the asymptotic limits of the Fcrmi- 
Dirac distribution (see Fig. 3): 

Non degenerate limit: In the limit A +00, the Lynden-Bell distribu- 
tion reduces to the MaxwcU-Boltzmann distribution 



/3 



1/2 



fQSsiP) - ( ^ ) e-^' (23) 
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C/2 




0.05 



Fig. 3. Spatially homogeneous Lyndcn-BcU distribution function for increasing values 
of A (top to bottom). For A = 0, the distribution reduces to a step function (completely 
degenerate) and for A — ► +oo, it becomes equivalent to the MaxwoU-Boltzmann distribu- 
tion (non degenerate). In the figure, we have taken /o = 0.13 and f3 has been calculated 
from Eq. (16). 



Since / <C /o, this corresponds to a dilute limit (or to a non degenerate 
limit if we use the terminology of quantum mechanics). The non degenerate 
limit corresponds, for a given value of /o, to A ^ +00, (3^0 and U — > +00. 
For /o +00, we are always in the non degenerate limit, for any /3 and U . 
In that case, the caloric curve takes the "classical" expression 

[/=— + -, (24) 
2/3 2 ^ ^ 

a relation that can be directly obtained from Eq. (23). For the water-bag 
initial condition, the non degenerate limit A +00 corresponds, for a 
given Mq, to U ^ +00. The non degenerate limit /o +00 corresponds 
to A9 ~ leading to Mq = 1 for any U In the non degenerate limit, the 



'^It also corresponds to Ap = leading to U = Umj{^{Mo) for any Mq. However, in that 
case the homogeneous phase does not exist since Umin{Mo) < Umin{Mo) so this case 
will not be considered here. 
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Lynden-Bell statistical equilibrium state describing the QSS has the same 
structure as the Boltzmann distribution describing the coUisional statisti- 
cal equilibrium state (but with, of course, a completely different interpre- 
tation). 

Completely degenerate limit: In the limit A ^ 0, the Lyndcn-BcU distri- 
bution (15) reduces to the Heavisidc function 

^'^SS{P) = ^, (25) 

where 

PF = (26) 
47r/o 

is a maximum velocity. The distribution (25) is similar to the Fermi distri- 
bution in quantum mechanics and pf is similar to the Fermi velocity. Thus, 
the limit A ^ corresponds to a completely degenerate limit in the quan- 
tum mechanics terminology. The completely degenerate limit corresponds 
to A — > 0, /3 -t-oo and U — ^ Uminifo) given by (21). This result can be 
directly obtained from Eq. (25). This is the minimum energy of the homo- 
geneous Lynden-Bell distribution for a fixed /q. This is also the minimum 
energy of any homogeneous distribution with / e {0, /□}. It corresponds 
to a water-bag initial condition with zero magnetization Mq = 0. If we 
start from a water bag initial condition, the completely degenerate limit 
corresponds to Mq = for any U and to U — Umini^Io) for any Mq- A 
stable water-bag initial condition with Mq = is a maximum Lyndcn-BcU 
entropy state, so it does not mix at all. 

In conclusion, we have in the general case, using the {U, /o) variables: 

• non degenerate limit 

— U +00 for any /o 

— fo +00 for any U 

• completely degenerate limit 

— U ^ Uminifo) for any fo 

For the water-bag initial condition, using the ([/, Mq) variables, we have: 

• non degenerate limit 

— ?7 ^ +00 for any Mq 

— Mq = 1 for any U 

• completely degenerate limit 
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- Mo = for any U 

- U = Umin{Mo) for any Mq 

4. Stability of the Lynden-Bell homogeneous phase 

We have seen that the maximization problem (6) provides a condition of 
thermodynamical stability (in Lyndcn-BcU's sense) and a condition of non- 
linear dynamical stability with respect to the Vlasov equation. We thus 
have to select the maximum of S at fixed U, A4. Indeed, a saddle point of 
S is unstable and cannot be obtained as a result of a violent relaxation. Let 
us consider the minimization problem 

m_in{i^[7] = t/[7] - r5[7] I M[f] = l}. (27) 
/ 

The criterion (6) can be viewed as a criterion of microcanonical stabil- 
ity and the criterion (27) as a criterion of canonical stability where F is 
interpreted as a free energy. Quite generally, the solutions of (27) are solu- 
tions of (6) but the reciprocal is wrong in case of ensemble inequivalence. 
Therefore, in the general case, criterion (27) forms a sufficient (but not 
necessary) condition of thermodynamical and formal nonlinear dynamical 
stability. In the present case, it can be shown that, if we restrict ourselves 
to spatially homogeneous solutions the ensembles are equivalent so that 
the set of solutions of (27) coincides with the set of solutions of (6) . There- 
fore, considering homogeneous states, criterion (27) forms a necessary and 
sufficient condition of thermodynamical and formal nonlinear dynamical 
stability. We shall therefore consider in the following the minimization prob- 
lem (27), which has been studied in^'^^ for general functionals of the form 
S[f] = — / C{f )d9dp where C is a convex function. A simple stability cri- 
terion has been obtained in the case where the steady state is spatially 
homogeneous, which can be expressed in terms of the distribution function 
as:22 

/■+°° f'(v) 

l + n ^-^dp > 0. (28) 

J-oo P 

It was shown in® that the same criterion can be expressed simply in terms 
of the density as 

cl>L (29) 



'^Concerning spatially inhomogeneous solutions, the microcanonical and canonical en- 
sembles are not equivalent in the region of first order phase transition. This important 
point will be further developed in a future contribution. 
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where = p'{p) is the velocity of sound in the corresponding barotropic 
gas. The equivalence between the criteria (28) and (29) is proven in.^ In 
this Section, we apply these criteria to the Lynden-Bell distribution (15). 
It is shown in^" that the criteria (28) and (29) can be rewritten: 



/_i/2(A)A|ri/2(A)|<^^-L_. (30) 



If the DF satisfies (30) <^ (28) ^ (29), then it is (i) Lynden-BcU thermody- 
namically stable (ii) formally nonlinearly dynamically stable. Otherwise, it 
is (i) Lynden-Bell thcrmodynamically unstable (ii) linearly dynamically un- 
stable.^ For given /o, the relation (30) determines the range of A for which 
the homogeneous distribution is stable/unstable. Then, using Eqs. (16), we 
can determine the range of corresponding energies. Specifically, the criti- 
cal curve Uc{fo) separating stable and unstable homogeneous Lynden-Bell 
distributions is given by the parametric equations 

/-V,(W1„.(A,| = a-i.^jM^ ,31) 

where A goes from (completely degenerate) to -l-oo (non degenerate). 
This leads to the phase diagram reported in Fig. 6. In fact, the criteria 
(28) and (29) only prove that / is a local entropy maximum at fixed mass 
and energy. If several local entropy maxima are found, we must compare 
their entropies to determine the stable state (global entropy maximum) and 
the metastable states (secondary entropy maxima). For systems with long- 
range interactions, metastable states have in general very long lifetimes, 
scaling like , so that they are stable in practice and must absolutely be 
taken into account. ^'^■^^ 

For /o —^ +O0, we are in the non degenerate limit A — > -foo and the 
stability criterion (30) for the homogeneous phase becomes 

f/ > t/c = J. (32) 

This returns the well-known nonlinear dynamical stability criterion (with 
respect to the Vlasov equation) of a homogeneous system with Maxwellian 
distribution function (sec, c.g.,^'^^). This also coincides with the ordinary 
thcrmodynamical stability criterion applying to the collisional regime, for 
t — > +00, where the statistical equilibrium state is the Boltzmann distribu- 
tion for / (without the bar!). 

On the curve U = U„mi{fo), we are in the completely degenerate limit 
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A ^ and the stability criterion (30) for the homogeneous phase becomes 

/o<(/o). = ^, ^.e U>U.^1-^. (33) 

This is the well-known nonlinear dynamical stability criterion (with respect 
to the Vlasov equation) of the water-bag distribution (see, e.g.,^'^^). 

If we start from a rectangular water-bag initial condition and use the 
(J7, Mq) variables, we must express /o in terms of U and 7\fo using Eq. (17). 
Then, the critical curve Uc{Mq) separating stable and unstable homoge- 
neous Lynden-Bell distributions is given by the parametric equations 



(34) 



" 48[(2t/-l)(A0)2+sin2A0]' ° M ' 

where A goes from (completely degenerate) to +oo (non degenerate). This 
leads to the phase diagram reported in Fig. 4. For Mq = 1, we get A +oo 
so we are in the non degenerate limit and the critical energy is Uc = 3/4. 
For Mo = 0, we get A = so we are in the completely degenerate limit and 
the critical energy is Uc = 7/12. 



5. The rectangular water bag initial condition: phase 
diagram in the (A^o, U) plane 

We first comment on the structure of the phase diagram in the {Mq,U) 
plane when we start from a water-bag initial condition. In Fig. 4 the tran- 
sition line Uc{Mo) divides the region of the plane where a homogeneous 
(-^QSS = 0) state is predicted to occur (upper area), from that (lower area) 
associated to a non- homogeneous phase (Afqss 7^ 0). Along the transition 
line two distinct segments can be isolated: the dashed line corresponds to 
a second order phase transition, the full line refers to a first order phase 
transition. First and second transition lines merge together at a tricritical 
point, approximately located at {Mo,U) = (0.17,0.61). The lateral edges 
of the metastability region associated to the first order transition line are 
also reported in the inset of Fig. 4. 

The correctness of the above analysis is assessed in^^ where numerical 
simulations of the HMF model (1) are performed for different values of 
the system size N. The transition predicted in the realm of Lynden Bell's 
theory are indeed numerically observed, thus confirming the adequacy of the 
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Homogeneous phase (stable) U^(M|j) ^^"^ 

Homogeneou s phase (unstable) 




'1 
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1 ' 














y ■< HI 

,1 


1 



0.2 



0.4 



M 



Forbidden region 







Fig. 4. Theoretical phase diagram in the control parameter plane {Mo,U) for a rect- 
angular water— bag initial profile. The dashed line Uc{Mo) stands for the second order 
phase transition, while its continuation as a full line refers to the first order phase transi- 
tion. The full dot is the tricritical point. The region below the thick full line Umin{!^o) 
corresponds to the forbidden domain of the parameter space. The region of existence of 
the homogeneous solution is delimited from below by the thin full line Umini^Io) (see 
Eq. (22)). Inset: zoom of the first order transition region. Dash-dotted lines represent 
the limits of the metastability region. In region (II), delimited from above by the upper 
dash-dotted line and from below by the full line, the homogeneous phase is fully stable 
and the inhomogeneous phase is metastable. In region (III), located below the full line 
and above the lower dash-dotted line, the homogeneous phase is metastable and the in- 
homogeneous phase is fully stable. These labels will appear again in Fig. 8, in connection 
with the discussion of the general case. 



proposed interpretative scenario. Moreover, the coexistence of homogeneous 
and inhomogeneous phases, corresponding to different local maxima of the 
entropy, is also reported in^^ where the probability distribution function of 
M is reconstructed. 

For all values of the energy larger than Uc{Mo) (where the homoge- 
neous phase is stable), the systems can potentially encounter a homoge- 
neous quasi-stationary phase which slows down the approach to the ther- 
modynamic equilibrium. Notice that above Uc = 0.75, the equilibrium value 
of the magnetization is also zero: there is hence no macroscopic transition 
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for M{t) of the type displayed in Fig. 1. One has therefore to rely on other, 
quantitative, indicators to monitor the dynamical state of the system, 
and eventually assess the presence of a QSS. This explains why in the past 
QSSs regimes where believed to be localized only in specific energy win- 
dows below the critical threshold Uc- Significant deviations from equilibrium 
are instead detected even for U > Uc, as reported in Fig. 5. Single parti- 
cle velocity distributions reconstructed from direct A^-body simulations at 
U ~ 0.85 display a bumpy profile, analogous to the one discussed in^^ for 
the reference energy U = 0.69 < Uc- Interestingly, for specific choices of 
the initial magnetization, the two bumps are even more pronounced than 
those analyzed in.^^ The presence of these bumps shows that relaxation 
is incomplete. These bumps correspond to the "vortices" (or phase space 
clumps) visible in Fig. 2. This state is stationary for the Vlasov equation, 
but Quasi-Stationary for the A^-body simulation. Hence, in the long run, 
the two "vortices" will merge, due to finite N effects, as the HMF system 
proceeds towards Boltzmann-Gibbs equilibrium. 



M =0.1 M =0.3 M =0.5 




Fig. 5. Probability distribution function of velocities f{p), for U=0.85 and different 
initial magnetization, as reported in the legend of each panel. 
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6. The general case: phase diagram in the {fQ,U) plane 

In the two-levels approximation, the Lynden-Bell equilibrium state depends 
only on two control parameters (U, fo). This is valid for any initial condi- 
tion with f{9,p, t = 0) G {0, fo}, whatever the number of patches and their 
shape. The variables (U, Mq) used in the previous section are valid only 
for a rectangular water-bag configuration. Furthermore, we note that two 
configurations with different values of (U, Mq) in a rectangular water-bag 
configuration can correspond to the same values of ([/, /o), hence to the 
same Lynden-Bell equilibrium state. To avoid any redundance, it is prefer- 
able to use the general variables {U, fo). Therefore, in the following, we shall 
discuss the general phase diagram in the (J7, fo) plane^*^ and compare it with 
the one in the {U, Mq) plane for the rectangular water-bag assumption. 

The stability diagram of the homogeneous Lynden-Bell distribution (15) 
is plotted in Fig. 6 in the {fo,U) plane. The representative curve Uc{fo) 
marks the separation between the stable (maximum entropy states) and 
the unstable (saddle points of entropy) regions. We have also plotted the 
minimum accessible energy for the homogeneous phase UmAnifo)- Here, the 
term "unstable" means that the homogeneous Lynden-Bell distribution is 
not a maximum entropy state, i.e. (i) it is not the most mixed state (ii) it is 
dynamically unstable with respect to the Vlasov equation. Hence, it should 
not be reached as a result of violent relaxation. One possibility is that the 
system converges to the spatially inhomogeneous Lynden-Bell distribution 
(15) with Mqss 7^ which is the maximum entropy state (most mixed) 
in that case. Another possibility, always to consider, is that the system 
docs not converge towards the maximum entropy state, i.e. the relaxation 
is incomplete?^ 

Let us enumerate some properties of the (/o, U) phase diagram of Fig. 6. 
For U > Uc = 3/4 (supercritical energies), the homogeneous phase is always 
stable (maximum entropy state) for any fo (recall that the curve Udfo) 
3/4 for fo +oo). On the other hand, there exists a critical point (we 
shall see later that it corresponds to the tricritical point of Fig. 4) located 
at 

(fo)* = 0.10947..., [/, = 0.608... (36) 

For fo < (/o)*, the homogeneous phase is always stable (maximum entropy 
state) for any U > C/,„i„(/o). For fo > (/o)c, the homogeneous phase is 
stable for U > Udfo) and unstable for Umin{fo) < U < Udfo)- This range 
of parameters corresponds to a second order phase transition. For [fo)* < 
fo < (/o)c, there is an interesting regime with a "re-entrant" phase^^ . 



February 2, 2008 12:57 



WSPC - Proceedings Trim Size: 9in x 6in marseillell 



19 



0.75 



0.7 



;d 0.65 - 



0.6 



0.55 



1 ' 


Homogeneous phase 
(stable) 




x»i - 

c 


. (fo)/ 


/ Homogeneous phase 
(unstable) 


- (fo)c 

1 


"^^^-^^ ^ (f^ 



0.1 



0.12 



0.14 



0.16 



Fig. 6. Stability diagram of the homogeneous phase in the (/o, (7) plane. The homoge- 
neous phase only exists above the line Uminifo)- It is stable above the line J7c(/o) a^nd 
unstable below it. For fixed fo G [(/o)*, {fo)c], there is a re-entrant phase as we progres- 
sively lower the energy: the homogeneous phase is stable for U > Uc (/o)> unstable for 
C/f'C/o) >U> ui^\fo) and stable again for r/<''(/o) >U> Urmnifo)- 



(2) (2) 

The homogeneous phase is stable for U > Uc {fo), unstable for Uc (/o) > 
U > J/i^H/o) a-nd stable again for Uc^\fo) > U > U,-nin{fo)- This range of 
parameters corresponds to a first order phase transition. 

To make the connection between the phase diagram (/o , U) obtained 
in^" and the phase diagram (Afg, U) obtained in,^^ we can plot the iso-Mo 
lines in the (/g, U) phase diagram. If we fix the initial magnetization Mq, 
or equivalently if we fix the parameter A0, the relation between the energy 
U and /o is 

1 1 /sinA6'\^ 1 



^^''^■^°)=6(4A0/o)2"2 +2- ^^^^ 



Therefore, the iso-Mo lines are of the form 



t/Ae(/o) = ^^-S(A^), (38) 
/o 
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Fig. 7. Iso-A/o lines in the {fo,U) phase diagram. This graphical construction allows 
one to make the connection between the (/o, U) phase diagram of Fig. 6 and the (Mq, U) 
phase diagram of Fig. 4. We can vary the energy at fixed initial magnetization by fol- 
lowing a dashed line. The intersection between the dashed line and the curve Uminifo) 
determines the minimum energy Umin{Mo) of the homogeneous phase. The intersection 
between the dashed line and the curve Uc{fo) determines the energy Uc{Mq) below which 
the homogeneous phase becomes unstable. 



with A{M) 



and B{M) 



1 Ae^ 



i which are easily rep- 



resented in the (/o, U) phase diagram (see Fig. 7). As an immediate conse- 
quence of this geometrical construction, we can recover the minimum energy 
of the homogeneous phase for a fixed initial magnetization A/q (or A9). In- 
deed, for a given A9, the homogeneous phase exists iff UAe{fo) > Umin{fo) 
leading to 



fo < (/o) 



487r2 sin^ AO ' 

This corresponds to ?7 > UminiMo) — C/Ae((/o)Ae) leading to 



rr rr r N ^ f Siu^ 



1 



(39) 



(40) 
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Fig. 8. Zones of mctastability in the (/o, U) phase diagram. In region (I), the homoge- 
neous phase is fully stable and the inhomogeneous phase is inexistent. In region (II) the 
homogeneous phase is fully stable and the inhomogeneous phase is metastable. In region 
(III) the homogeneous phase is metastable and the inhomogeneous phase is fully stable. 
In region (IV) the homogeneous phase is unstable and the inhomogeneous phase is fully 
stable. The three curves separating these regions connect themselves at the tricritical 
point. In region (V) the homogeneous phase is inexistent. The corresponding caloric 
curves as well as the absolute minimum energy UminHo) of the inhomogeneous phase 
will be determined in a future contribution. 

which is identical to (22). Note again that A6 is related to A/q by Eq. (14). 
Figure 7 is in good agreement with the structure of the phase diagram 
in the {U, Mq) plane. Indeed, along an iso-Mg line, we find that for large 
energies U > Uc{Mo) the homogeneous phase is stable and for low energies 
U < Uc{Mo) the homogeneous phase becomes unstable. In that case, there 
is no re-entrant phase. 

In,^" only the stability of the homogeneous phase has been studied, i.e. 
whether it is an entropy maximum at fixed mass and energy or not. The 
question of its mctastability, i.e. whether it is a local entropy maximum with 
respect to the inhomogeneous phase, has not been considered. However, 
considering Fig. 7 and comparing with the results of,^^ we conclude that 
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there must exist zones of metastability in the (/o, U) phase diagram. They 
have been represented in Fig. 8. In region (I) the inhomogeneous phase 
does not exist while the homogeneous phase is fully stable. In region (II), 
the inhomogeneous phase appears but is mctastablc while the homogeneous 
phase is fully stable. In region (III), the inhomogeneous phase becomes fully 
stable while the homogeneous phase becomes metastable. In region (IV), 
the homogeneous phase becomes unstable while the inhomogeneous phase 
is fully stable. The three curves separating these regions connect themselves 
at a tricritical point. This is clearly the same as in Fig. 4. Using Eqs. (36), 
(37), (14) we find that it corresponds to 

[/* = 0.608..., (A/o), = 0.1757... (41) 

with = 2.656.... Therefore, the phase diagrams in (/o, U) and {Mq, U) 
planes are fully consistent. Note, however, that the physics is different 
whether we vary the energy at fixed /o or at fixed Mq. In particular, there 
is no "re-entrant" phase when we vary the energy at fixed Mq^^ while a 
"re-entrant" phase appears when we vary the energy at fixed /o.^° 

7. Conclusions 

In this paper, we have discussed the emergence of out-of-equilibrium Quasi 
Stationary States (QSSs) in the Hamiltonian Mean Field (HMF) model, a 
paradigmatic representative of systems with long-range interactions. The 
analysis refers to a special class of initial conditions in which particles are 
uniformly occupying a finite portion of phase space and the distribution 
function takes only two values, respectively and /q. The energy can be 
independently fixed to the value U . 

The Lynden-Bell maximum entropy principle is here reviewed and 
shown to result in a rich out-of~equilibrium phase diagram, which is conve- 
niently depicted in the reference plane (/o, U)}^ When considering a rect- 
angular water-bag distribution the concept of initial magnetization, Mq, 
naturally arises as a control parameter and the different QSSs phases can 
be represented in the alternative space (Mq, U)}"^ In both settings first and 
second order phase transitions are found, which merge together in a tri- 
critical point. These findings have been tested versus numerical simulation 
in,^^ where the adequacy of Lynden-Bell theory was confirmed. 

A formal correspondence between the two above scenarios is here drawn 
and their equivalence discussed. It is worth mentioning that swapping from 
one parametric representation to the other allows us to put the focus on in- 
triguingly different physical mechanisms, as it is the case of the "re-entrant" 
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phases discussed in Section 6. 

Further, in this paper, we have provided an analytical characterization of 
the domain of existence of the Lynden-Bell spatially homogenous phase and 
investigated its stability. Homogeneous QSS are also expected to occur for 
U > Uc = 3/4, a claim here supported by dedicated numerical simulations. 

Despite the fact that Lyndcn-BcU's theory results in an accurate tool to 
explain the peculiar traits of QSSs in HMF dynamics, one should be aware 
of the limitations which are intrinsic to this approach. Most importantly, 
Lynden-Bell's recipe assumes that the system mixes well so that the hy- 
pothesis of ergodicity, which motivates the statistical theory (maximization 
of the entropy), applies. Unfortunately, this is not true in general. Several 
example of incomplete violent relaxation have been identified in stellar dy- 
namics and 2D turbulence (see some references in^^) for which the QSSs 
cannot be exactly described in term of a Lynden-Bell distribution. Also in 
this case, however, the QSSs are stable stationary solution of the Vlasov 
equation and novel analytical strategies are to be eventually devised which 
make contact with the underlying Vlasov framework. 
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